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Abstract 

\ 

In  this  report,  we  prove  that  the  state  vector  of  a  two-point  boundary 
value  model  is  a  reciprocal  process. 
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1.  Introduction 


When  modeling  physical  systems  in  which  one  or  more  of  the 
independent  variables  is  spatial,  the  a  priori  information  is  usually  given  in 
the  form  of  boundary  conditions  rather  than  initial  conditions.  Such  systems 
arise  in  acoustics  and  oceanography,  for  example.  Consider  the  following 
linear,  two-point  boundary  value  model  for  ^  <  t  <  tx : 


x{t)  =  Ax[t)  +  Bu[t) 

(1) 

r  =  V0x[t0)  +  V1x[t1) 

12) 

where  u  is  a  Gaussian,  zero-mean,  unit  intensity  white  noise  with  m 
components,  r  is  a  Gaussian,  zero-mean  random  vector  with  n  components 
and  is  independent  of  u,  x  is  the  n-component  state  vector,  and  A,  B,  V0,  V1 
are  appropriate  matrices. 

Structural  properties  of  such  two-point  boundary  value  models  have 
been  studied  in  [5]-[6].  Estimation  problems  based  on  such  models  have 
been  treated  in  [2]-[3],  [7]-[8],  As  is  customary,  we  will  assume  that  the 
model  is  well-posed ;  i.e.,  that  u  and  r  give  rise  to  a  unique  x.  This  will  be 
the  case  if  and  only  if  the  matrix  V0  +  Vj<2>(ti  ,£<,)  is  nonsingular,  in  which  case 
we  can  assume,  without  loss  of  generality,  that 


Vo+V^MoW 

where  0  is  the  state  transition  matrix  of  (1). 

Since,  in  general,  xfto)  is  correlated  with  u,  the  state  vector  x  is  not 
Markovian.  However,  as  we  shall  prove,  it  does  possess  the  weaker  property 
of  reciprocity,  which  will  be  defined  presently.  First,  we  note  that  Krener 
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stated  this  result  in  [5]  and  sketched  a  method  of  proof  for  it.  However,  we 
have  been  unable  to  complete  the  proof  along  the  suggested  lines,  and  were 
thus  motivated  to  find  a  different  approach. 

A  process  is  reciprocal  on  an  interval  if,  given  any  subinterval,  the 

values  of  the  process  inside  the  subinterval  are  conditionally  independent  of 

the  values  outside  the  subinterval,  given  the  values  of  the  process  at  the  two 

endpoints  of  the  subinterval.  For  information  on  reciprocal  processes,  see 

[1],  [4]-[5]  and  the  references  therein.  In  a  recent  paper  [4],  Carmichael  et 

al  give  a  very  useful  characterization  of  reciprocal  processes,  based  on  the 

projection  theorem.  They  prove  that  a  process  x  is  reciprocal  if  and  only  if, 

for  all  ae  [io  ,  tj,  the  two  (one-sided)  error  processes  e1  and  e,  given  by 

ex  [t)  =  x(t)  -  E[x(t)lx(a)]  for  tQ  <  t  <  a 
e2{t)  =  x[t)~ E[x(f )l x(a)]  for  a  <t<  tj 

are  Markovian.  Our  proof  is  based  on  this  result,  and  on  the  fact  that  if  a 
Gaussian  process  has  a  semi-separable  covariance  function,  it  is  Markovian. 


2.  The  Main  Result 


The  boundary  value  model  (1),  (2)  can  be  solved  for  x(a),  for  all 
as  [<o  ,  tj,  as  follows: 

x(a)  =  0(a,to)r+  ftlG(a,T)Bu(T)dr 

■'to 

where  G  is  the  Green's  function  given  by 

f-0(a,to)V1<J>(t1,t)  for  a  <  t 


G(a,t)  = 


<Z>(a,t0)V0<J>(t0,t)  for  a  >  t 


(3) 
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Let  R[a,a)  denote  the  covariance  matrix  of  xfa).  We  assume  this  matrix  is 


nonsingular  for  all  ae  ,  tj.  From  (1)  and  the  definitions  of  ex  and  e2  given 
above,  it  is  clear  that  for  4  <  t  <  a 


e1[t)  =  Ae1{t)  +  Bu{t),  e1[a)  =  0 

(4) 

and  for  a  <  t  <  tx 

e2(t)  =  Ae2 (t)  +  Bu(t),  e2 (a)  =  0 

(5) 

where 

u(t)  =  u(t)-£[u(t)lx(a)] 

=  u(t)  -  E[i4t)x'(a)]£" 1  (a,  a)x(a) 

=  u[t)  -  B'G'{a,t)R~Ha,a)x{a) 

and 

E[u(£)u'(s)]  =  I8[t  -  s)  -  BTT(a,£)£-1(a,a)G(a,s)B 

(6) 

Theorem: 

The  state  vector  process  x  is  reciprocal. 

Proof:  As  discussed  above,  it  is  sufficient  to  prove  that  e1  and  e2 
are  Markovian  for  all  a  e  ,  tj.  Consider  first  the  case  of  e2.  For  a<t<tx 
and  a  <  s  <  tx ,  its  covariance  function  is  given  by 
Kz(t,s)  =  E[e2[t)e'2{s)] 

=  £  £<P(f,  r)B{E[u[r)u'{(7)]}B'0'{s,  o)dadx 
where  we  have  used  (5).  Substituting  in  the  above  using  (6)  produces 
KJt,s)  =  f  (P{t,T)BB'<P'(s,T)dx 

Ja 

-[£  <P[t,  x)BB'G'{a,  t)  drjf?-1  (a,  a)[j*G(a.  t)BB'<Z>'(s,  t)  dr 
and  using  equation  (3)  for  the  Green's  function,  we  get 
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K2(t,s)  =  [min<t‘s>  0(t,T)BB'0'(sfT)dr 

•  a 

-[£  0{t, z)BB'®%,z)dz  V{0'(a.to)R-l(a,a)<P[a,to)V1  [ £<£(£, , x)BB'&'{s,  x)dx 

f C2(t)D2{s )  for  a  <s<t<tj 
[D2  ( t)C2  (s)  for  a  £  t  <  s  <  tj 

where 

C2(f)  =  ^(f.t1)-[£0(t,T)BB,<2>/(t1,r)dT]v/<J>/(a,to)J?-1(a,a)^(a,to)V1 
D2  (t)  =  £  <Z>(tj ,  r)BB'0'[t,  f)  dr 

Since  its  covariance  function  is  semi-separable,  e2(t),  a<t<tlt  is  Markovian 
for  all  as  ,  t,].  Using  the  same  approach,  it  can  be  shown  that,  for 
*o  <  t  <  a  and  ^  <  s  <  a, 

K1[t,s)  =  E[el{t)ei{s)] 

=  fa  0{t,t)BB'0'{s,r)dz 

Jma x(t.s) 

-[£<P(t,  z)BB'&'(t0,  t)dry'0'(ado)R-Ha,a)0{ado)Vo^\y[to,  z)BB’<I>'{s,  x)dx 

_(C1{t)D1{s)  for  t0<t<s<a 
[Di(t)Ci(s)  for  t0  S  s  <  t  <  a 

where 

CJt)  =  ^>(t,t0)  -  ^(t.rJBB'^'Uo,  T)dTjvo'0'(a,to)K  ■1(a,a)‘3>(a,t0)V0 
D1(t)  =  j0{to,x)BB'0'{t,x)dz 

Therefore,  ejt),  t<  a,  is  Markovian  for  all  ae  [^  ,  tj,  and  the  theorem  is 
proved. 


5 


References 


[1]  J.  Abrahams  and  J.  B.  Thomas,  "Some  comments  on  conditionally 
Markov  and  reciprocal  Gaussian  processes,"  IEEE  Transactions  on 
Information  Theory,  vol.  27,  pp.  523-525,  1981. 

[2]  M.  B.  Adams,  A.  S.  Willsky,  and  B.  C.  Levy,  "Linear  estimation  of  boundary 
value  stochastic  processes,"  IEEE  Transactions  on  Automatic  Control,  vol. 
29,  pp.  803-821,  1984. 

[3]  A.  Bagchi  and  H.  Westdijk,  "Smoothing  and  likelihood  ratio  for  Gaussian 
boundary  value  processes,"  IEEE  Transactions  on  Automatic  Control,  vol.  34, 
pp.  954-962,  1989. 

[4]  J.-P.  Carmichael,  J.-C.  Masse,  and  R.  Theodorescu,  "Representations  for 
multivariate  reciprocal  Gaussian  processes,"  IEEE  Transactions  on 
Information  Theory,  vol.  34,  pp.  155-157,  1988. 

[5]  A.  J.  Krener,  "Reciprocal  processes  and  the  stochastic  realization 
problem  for  acausal  systems,"  Modelling,  Identification,  and  Robust  Control, 
C.  I.  Byrnes  and  A.  Lindquist,  Eds.,  Amsterdam:  Elsevier,  1986,  pp.  197-209. 

[6]  A.  J.  Krener,  "Acausal  realization  theory  I:  Linear  deterministic  systems," 
SIAM  J.  Control  and  Optimization,  vol.  25,  pp.  499-525,  1987. 

[7]  L.  R.  Riddle  and  H.  L.  Weinert,  "Fast  algorithms  for  the  reconstruction  of 
images  from  hyperbolic  systems,"  Proc.  IEEE  Conference  on  Decision  and 
Control,  1986,  pp.  173-178. 

[8]  L.  R.  Riddle  and  H.  L.  Weinert,  "Recursive  linear  smoothing  for 
dissipative  hyperbolic  systems,"  Mechanical  Systems  and  Signal  Processing, 
vol.  2,  pp.  77-96,  1988. 


6 


ia  report  security  classification 

Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


ID.  RESTRICTIVE  MARKINGS 


2D.  OECLASSIFICATION/DOWNGRAOING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBERIS) 

JHU/ECE  89-14 


6a.  NAME  OF  PERFORMING  ORGANIZATION  Sb.  OFFICE  SYMBOL 

(If  applicable  / 

The  Johns  Hopkins  University 


6c.  AOORESS  (City.  State  and  ZIP  Code) 

Baltimore,  Maryland  21218 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBERIS) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Office  of  Naval  Research 


7b.  AOORESS  (City.  State  and  ZIP  Code) 

Arlington,  Virginia  22217 


a*.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 


Be.  AOORESS  (City.  State  and  ZIP  Code) 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

N00014-89-J-1642 


10.  SOURCE  OF  FUNDING  NOS 


PROJECT 

TASK 

WORK  UNIT 

NO. 

NO. 

NO. 

12.  PERSONAL  AUTHOR(S) 

Howard  L.  Weinert 


13a  TYPE  OF  REPORT 

Technica 


19.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVEREO 
FROM _  TO 


14.  OATE  OF  REPORT  IYr..  Mo..  Day)  15.  PAGE  COUNT 

1989,  September  22  6 


COSATI  COOES 


FIELO  GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverie  if  ne canary  and  identify  by  block  number) 

^  Boundary  value  models,"  reciprocal  processes  ,  ./  } 


19.  ABSTRACT  (Continue  on  reverie  if  neceuary  and  identify  by  block  number/ 

^  .  _ _  l 

In  this  report^  we  prove/ that  the  state  vector  of  a  two-point  boundary  value  model  is 
a  reciprocal  process. /  V',  ,,  ,  ,  y.r 


20.  OISTRIBUTION/AVAILABILITY  of  abstract 

unclassifieo/unlimiteo  S  same  as  rpt.  □  otic  users  G 


22a.  name  of  responsible  inoivioual 

Dr .  Julia  Abrahams 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22b.  TELEPHONE  NUMBER  22c.  OFFICE  SYMBOL 

(Include  Area  Code / 

(202)  696-4320 


DO  FORM  1473,  83  APR 


EDITION  OF  1  JAN  73  IS  OBSOLETE. 


